Mapping the local density of optical states of a photonic crystal with single quantum 

dots 



O 

(N 



Oh; 

a ■ 



> 

On 

o 



X 



Qin Wang,^ S0ren Stobbe,^ and Peter Lodahl-'^ 

^DTU Fotonik, Department of Photonics Engineering, 
Technical University of Denmark, 0rsteds Plads 343, 2800 Kgs. Lyngby, Denmark. 

We use single self-assembled InGaAs quantum dots as internal probes to map the local density 
of optical states of photonic crystal membranes. The employed technique separates contributions 
from non-radiative recombination and spin-flip processes by properly accounting for the role of the 
exciton fine structure. We observe inhibition factors as high as 55 and compare our results to local 
density of optical states calculations available from the literature, thereby establishing a quantitative 
understanding of photon emission in photonic crystal membranes. 

PACS number(s): 42.50.Ct, 78.67.Hc, 78.47.D- 



Photonic crystals (PCs) are artificial periodic dielec- 
tric materials that were originally proposed as a way 
to dynamically control spontaneous emission of light po- 
tentially leading to efficient light sources and solar cells 
Embedding light sources inside photonic crystals 
has enabled the experimental demonstration of modified 
spontaneous emission using dye molecules quan- 
tum wells or quantum dots [^-Q. The complexity 
of light-matter interaction in PCs is apparent since the 
experiments combine inhomogeneous dielectric materials 
varying on the scale of the wavelength with inherently 
mesoscopic quantum emitters. Within the validity of 
the dipole approximation, which is valid for standard- 
sized quantum dots in dielectric structures , the light- 
matter coupling is determined by two quantities: i) the 
transition dipole moment of the optical transition and 
ii) the local density of optical states (LDOS) projected 
onto the orientation of the dipole [ll|, The latter 

is a property of the PC that accounts for all available 
Bloch modes at a specific position and has been exceed- 
ingly challenging to calculate in realistic structures rele- 



vant for experiments ll|, |13[ . Consequently, a thorough 



understanding of the LDOS and therefore the potential 
of spontaneous emission control has been lacking in PCs. 
Here we present the experimental mapping of the LDOS 
of PCs by time-resolved emission studies of single quan- 
tum dots (QDs) with well-characterized optical proper- 
ties. 

The focus of the present work is to use QDs as inter- 
nal probes of the LDOS. This quest requires a detailed 
understanding of the optical properties of QDs in inho- 
mogeneous media and in particular the role of the exciton 
fine structure. Our previous work on dielectric interfaces 
demonstrated the necessity to consider the role of dark 
excitons and proved the validity of the dipole approx- 
imation for standard-sized QDs in dielectric structures 
[loj . while a breakdown of the dipole approximation was 



observed near metallic interfaces (15[. Building on this 
knowledge, we present a new and general method to ex- 
tract the radiative decay rate of single QDs as opposed to 
the total decay rate that is otherwise directly measured 



in time-resolved spectroscopy and is significantly infiu- 
enced by dark excitons and other non-radiative decay 
processes. From the radiative decay rate we determine 
the LDOS and experimentally map it by recording decay 
curves of many single QDs positioned throughout the 2D 
PC membranes. We compare our results with LDOS sim- 
ulations performed with 3D finite-difference time-domain 
(FDTD) simulations available in the literature [l^ and 
explicitly demonstrate the importance of extracting the 
radiative decay rate. We observe record high sponta- 
neous emission inhibition factors of 55 compared to QDs 
in homogenous media, which proves the potential of 2D 
PC membranes for applications where spontaneous emis- 
sion is a nuisance. The frequency dependency of the 
LDOS is mapped out both inside and outside the 2D 
band gap, and in the latter case also enhancement of the 
spontaneous emission rate is observed. 

Within the dipole approximation and for weak light- 
matter interaction strengths, the radiative decay rate 
is directly proportional to the LDOS: 7iad(ro,w) — 
|/x|^p^(ro,w), where p;^(ro,a;) is the projected 
LDOS evaluated at the position ro and emission fre- 
quency uj of the emitter, eq is the vacuum permittiv- 
ity, and fi is the transition dipole moment. The LDOS 
describes the electromagnetic environment and can be 



calculated from the dyadic Green's function [16|. Suc- 
cessfully extracting the radiative decay rate provides a 
mean to obtaining the LDOS by comparing to the radia- 
tive decay rate in a homogeneous medium '^]^^^{ti') using 
the relation 



Pp(i"o,w) 



7rad(ro,ll^) 



(1) 



where 



is the projected LDOS for a homoge- 



neous medium with refractive index n and c is the speed 
of light in vacuum. 

A detailed understanding of the exciton fine structure 
and the role of non-radiative recombination is required 
in order to use QDs as LDOS probes 0,113, El- An 
InGaAs QD has two optically active bright states with 
total angular momentum Jz = ±1 and two dark states 



with Jz = ±2, where the quantization axis z is the growth 
direction [001] [l8|. Due to the reduced symmetry and 
anisotropic exchange interactions, the two bright states 
are separated in energy and form two eigenstates X and 
Y named according to their dipole orientations ([110] or 

[110]), where \X\ = (|+1) + |-1))/V2, \Y\ = (|+1) - 
| — 1))/\/2. Similarly, the two dark states are separated 
into \X)^ = (1+2) + |-2))/V2, |r)d = (1+2) - \-2))/V2. 
The PC patterns can be aligned relative to X and Y, as 
schematically indicated in Fig. 1(a). 

The two bright and two dark states together with the 
ground state form a five-level model. \X)^^ couples to 
\X)^ through either an electron or a hole spin flip me- 
diated by phonons and exchange interaction [l^. We 
note that spin-flip processes coupling bright excitons, i.e. 
\X)^ ^ |F),^ are slow compared to the other decay pro- 
cesses and therefore can be abandoned in the analysis, as 
theoretically predicted and experimentally con- 

firmed from the large anisotropy in the decay rate for 
X and Y states observed in a PC jllj. Consequently, 
the five-level scheme can be simplified to the three-level 
system indicated in Fig. 1(c). The bright state can de- 
cay either through radiative or non-radiative processes 
with rates jf^'^ and T^lJd' respectively, and is coupled 
to the dark state through a spin-flip rate 7j^. Radia- 
tive transitions from the dark state to the ground state 
are forbidden, but non-radiative recombination is possi- 
ble with the rate 7j^a'd together with spin flips to the 
bright state with a rate 7^. Solving the resulting rate 
equations we find that the population of bright exci- 
tons decay bi-exponentially {t) — Afe^'^f * + Ase~^<' *, 

7fb + \/(7ii)' + (7id /2)^ 

- 
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Similar expressions 



can be obtained for p^{t). (0) and p^{0) are the ini- 
tial populations of the bright and dark states. For non- 
resonant and weak pumping we have p^ (0) = Pb"(0) = 
0.5, where we note that the presented results of the LDOS 
mapping are robust to deviations from this equilibrium 
condition. Furthermore, we have assumed: (i) 7jJ^ = 7^ 
and (ii) 7,'J^ad — '^mtd- (0 ^ good approximation for in- 
termediate temperatures (T = 10 K in our experiment), 
because 7bd — e^^'^^'^^'^^dh, and ksT ^ Shd, where fee is 
the Boltzmann constant, and dhd is the energy splitting 
between the bright and dark states (typically, a few hun- 
dred peV llSj). Assumption (ii) has been proven valid 
by experiments on QDs in dielectric media with a known 
LDOS lj,ll7|. The intensity measured in a time-resolved 
experiment is (t) — C'cT^d /^b' '^^tiere Cq is pro- 
portional to the total collection efficiency of the exper- 
imental setup, which is an overall scaling factor of the 
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FIG. 1. Fig. 1 (Color online) (a) Sketch of a PC membrane, 
top view (left) and side view (right). The orientations of the 
X and Y dipoles together with the sample growth direction 
(Z) are schematically indicated. The inset shows the first 
Brillouin zone with indications of the relevant symmetry di- 
rections; (b) A photonic-band diagram calculated for optical 
modes with in-plane polarization, where the gray region rep- 
resents the 2D photonic band gap (PBG) and "leaky" refers 
to the regions where light is not confined to the membrane; (c) 
Measured decay curves of QD A and B that are, respectively, 
tuned in- or outside the band gap, as schematically indicated 
in (b). The insets show the three- level diagram of the QD 
(left) and a measured emission spectrum (right). 



decay curves. After fitting each decay curve with a bi- 
exponential function, we obtain the four parameters 7f , 
-fl, A\, and Ai (for i = X,Y) and 7*^^, 7^,.^^ and 7^^ 
can be extracted from the relations described above. The 
amplitudes need to be corrected to account for residual 
population of the QD when it is re-excited by a light pulse 
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141: A\ = ^J[l-exp(-7fV)] and = ^][l-exp(-7|r)], 
where r is the excitation period, and A\ {A\) refers to 
the measured amplitude. We stress that the presented 
method is completely general, i.e., no assumptions about 
the magnitude of the rates have been implemented, and 
therefore InGaAs QDs can be employed as LDOS probes 
in any nanophotonic environment. 

We have carried out time-resolved measurements on 
single self-assembled QDs positioned in- or outside 2D 
PC membranes. The experiments are done in a flow 
cryostat at 10 K and the sample consists of a series of 
2D GaAs (n=3.5) PC membranes with a layer of InGaAs 
QDs (density ~ 80 /J,m~^) embedded in the center. The 
dimension of each PC membrane is 12 x 12 /xm^ with a 
thickness of 154 nm as shown in Fig. 1(a) where also the 
orientations of X and Y dipoles are indicated. Fig. 1(b) 
shows the photonic band diagram of the structure with 
indications of the 2D photonic band gap region and the 
continuum due to coupling to leaky modes out of the 
membrane. The lattice parameter (a) ranges from 200 to 
385 nm in steps of 5 nm, and the r /a ratio is fixed at 0.30, 
where r is the radius of the air hole. The QDs are excited 
with a PicoQuant PDL-800 pulsed diode laser at 781 nm 
with varying repetition rate (5, 10, 20, 40 MHz). Under 
a weak excitation condition, we identify neutral excitons 
by their excitation power and polarization dependence, 
i.e., multi-exciton complexes are excluded due to the ob- 
served linear power dependence before saturation, while 
single-charge excitons are found to be mono-exponential 
and have a very weak polarization dependence. We se- 
lect only QDs that emit within a narrow spectral range 
of 970 ± 5 nm, in order to probe QDs with similar oscilla- 
tor strengths [13] ■ In total 88 QDs in the PC are probed 
in addition to 5 QDs positioned outside the PC pattern 
that serve as reference. For each QD, a polarizer is used 
to record emission from either the X or the Y exciton 
state. 

Figure 1(c) shows two typical decay curves of QD A 
and B (X states) that are tuned respectively in- or out- 
side the band gap. Clearly, the decay of QD A is strongly 
suppressed as a consequence of the strongly suppressed 
LDOS associated with the 2D band gap of the structure. 
QD B is positioned in a PC that is designed to probe 
the LDOS at the edge of the band gap, and indeed a 
much faster decay is observed in this case. Based on the 
method explained above, we can extract the radiative de- 
cay rate of each exciton state and obtain the projected 
LDOS. The results are shown in Fig. 2 including a com- 
parison to the theory of Koenderink et al. [3] who cal- 
culated the projected LDOS for 4 specific positions in a 
PC membrane. In order to compensate for a slight dif- 
ference in membrane thickness between experiment and 
theory, we have extended the band gap width according 
to the theory of Ref. [22] ■ We observe a pronounced 2D 
band gap leading to a wide range of frequencies where the 
LDOS is strongly suppressed for both X and Y exciton 
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FIG. 2. Fig. 2 (Color online) Measured projected LDOS of 
PC membranes versus normalized emission frequencies for X 
and Y dipole orientations (data points). The dashed lines 
represent the calculated projected LDOS as taken from Ref. 
[13 ] . For reference, also the DOS of homogeneous GaAs is in- 
dicated (black dash-dotted lines) and the gray areas represent 
the band gap. 



states. Large point-to-point fluctuations in the experi- 
mental data are observed since the QDs exhibit different 
radiative decay rates due to their varying positions in the 
PC membrane, thus reflecting the sensitive spatial vari- 
ation of the LDOS. Outside the band gap an enhanced 
LDOS is observed relative to the DOS of a homogeneous 
medium. Consequently, the decay rate of QDs tuned to 
the band edge can be Purcell enhanced in this case by 
coupling to extended Bloch modes as opposed to local- 
ized cavity modes, which is the most common way of 
realizing the Purcell effect. In general good agreement 
between experiment and theory is observed, especially 
taking into consideration that the spatial sampling used 
in the theory is rather sparse, while in the experiment 
the QDs occupy all different positions in the high refrac- 
tive index material. Interestingly, the experimental data 
appear to be biased such that the measured LDOS seems 
systematically larger than predicted by theory, in partic- 
ular for normalized frequencies below the 2D band gap. 
This observation is likely an effect of unavoidable fabri- 
cation disorder that would be more severe for the small 
feature sizes corresponding to reduced frequencies below 
the band gap. Disorder has been found to lead to signif- 
icant modiflcations of the LDOS in PC waveguides 
and it is therefore important to quantify the role of dis- 
order in any PC application, which is done here through 
the comparison with theory. 

For reference, we have also measured the radiative de- 
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FIG. 3. Fig. 3 (Color online) Inhibition factors for five 
strongly inhibited QDs due to the 2D photonic band gap. 
The blue circular points are obtained from the full model ac- 
counting for the exciton fine structure and non-radiative re- 
combination, while the red diamond points are obtained by 
incorrectly using directly measured total decay rates. The 
dashed blue line shows the maximum values from simulations 

M. 



cay rate, non-radiative decay rate and spin-flip rate of 5 
QDs positioned outside the PC area and obtain the av- 
erage values 7^^^ = 1-1 ± 0-1 ns-\ 7^^^^ = 0-06 ± 0.05 
ns~"^ and 7^^ = 0.005 ± 0.002 ns~^ respectively, vifhich 
agree well with previous results l3, 14, 17|. Focusing on 
QDs positioned in the band gap we extract the inhibition 
factor relative to the reference value in the homogenous 
medium, see Fig. 3. A maximum inhibition factor of 55 is 
observed, which is the highest value ever reported in any 
PC. We stress the necessity of employing the presented 
method that accounts for the QD fine structure and non- 
radiative recombination in order to correctly extract the 
inhibition factors. Thus, Fig. 3 also displays the inhibi- 
tion factors derived by using the common but incorrect 
assumption that the directly measured total decay rate is 
dominated by radiative recombination. This comparison 
clearly illustrates the importance of employing the cor- 
rect microscopic model of the quantum emitter in order 
to use them as LDOS probes. 

In conclusion, we have presented a method to probe 
the LDOS of any nano environment by employing self- 
assembled InCaAs QDs. By properly accounting for the 
exciton fine structure, it is possible to extract the ra- 
diative decay rate and therefore eliminate effects from 
non-radiative recombination and spin-fiip processes. We 
presented a detailed frequency map of the LDOS and a 
detailed comparison to existing theory. Inhibition factors 
as high as 55 were observed inside the 2D band gap, thus 
clearly demonstrating the potential of PC membranes for 
efficient spontaneous emission inhibition. Our work is 



expected to lay the foundation for further exploitations 
of photonic crystal membranes for all-solid-state quan- 
tum electrodynamics experiments, where the LDOS is 
the essential quantity that controls not only spontaneous 
emission, but also, e.g., the Lamb shift 2J] or Casimir 
forces. 
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